Electron transport through interacting quantum dots 
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We present a detailed theoretical investigation of the effect of Coulomb interactions on electron 
transport through quantum dots and double barrier structures connected to a voltage source via 
an arbitrary linear impedance. Combining real time path integral techniques with the scattering 
matrix approach we derive the effective action and evaluate the current-voltage characteristics of 
quantum dots at sufficiently large conductances. Our analysis reveals a reach variety of different 
regimes which we specify in details for the case of chaotic quantum dots. At sufficiently low energies 
the interaction correction to the current depends logarithmically on temperature and voltage. We 
identify two different logarithmic regimes with the crossover between them occurring at energies of 
order of the inverse dwell time of electrons in the dot. We also analyze the frequency-dependent shot 
noise in chaotic quantum dots and elucidate its direct relation to interaction effects in mesoscopic 
electron transport. 



I. INTRODUCTION 

Low temperature electron transport in disordered con- 
ductors with electron-electron interactions remains one of 
the most intriguing topics in modern mesoscopic physics. 
Interplay between charge discreteness, quantum coher- 
ence, scattering and interactions yields a rich variety 
of non-trivial effects, many of which can adequately be 
described only within a rather complicated theoretical 
framework. 

One of the most successful theoretical approaches in 
mesoscopic physics is the scattering matrix formalism-jS.. 
In the absence of interactions this method allows for a 
complete and physically transparent analysis of electron 
dynamics in coherent conductors which includes not only 
electrical conductance but also shot noise'^ and eventu- 
ally all higher cumulants of the current operator—. For- 
mally it is possible to generalize the Landauer formal- 
ism to systems with interactionsiSi. Within this approach, 
however, the electron Green functions in the interacting 
region still remain to be evaluated diagrammatically or 
by other means. 

An alternative way to tackle the problems with dis- 
order and interactions is to employ non-perturbative 
path-integral-based techniquesSii. Recently it was 
demonstrated2i2ii£ that the advantages of the path in- 
tegral methods and the scattering matrix approach can 
be conveniently combined within one formalism, thereby 
providing a powerful tool to theoretically describe inter- 
action effects in mesoscopic conductors. Several impor- 
tant results have already been obtained in this way. For 
instance, with the aid of the instanton technique it was 
shown^ that - similarly to tunnel junctionsii*i^ - arbi- 
trary coherent conductors can exhibit the phenomenon 
of weak charge quantization. Unfortunately, in the limit 
of large dimensionless conductances g ^ I this effect is 
rather weak and it gains importance only at exponen- 
tially small temperatures and voltages. Furthermore, 
weak charge quantization vanishes completely even at 
T = provided at least one of the conducting channels is 



fully transparentSii^. The results^ were later confirmed 
and extended in Refs. ^3 andlisl 

Another - much more robust - effect of electron- 
electron interactions is the so-called interaction correc- 
tion to the / — U curve. This interaction correction 
is negative and its magnitude scales linearly with the 
parameter— 



(1) 



where T„ is the transmission of the n-th conducting chan- 
nel. In contrast to the effect of weak charge quantization, 
even for large conductances 5 3> 1 the interaction correc- 
tion remains clearly observable up to high temperatures 
and it vanishes only for P —^ 0, i.e. provided all the 
conducting channels in the system are fully transparent. 
Moreover, as both temperature and voltage get lowered, 
the interaction correction grows logarithmically-Ssi^ . As 
a result, at sufficiently low energies the smallness ^ 1/g 
gets compensated by a large logarithm and the system 
enters a non-perturbative regime where terms of all or- 
ders in the interaction need to be evaluated. Let us also 
recall that for a particular case of highly conducting tun- 
nel junctions this logarithmic enhancement of the inter- 
action correction is already well known for single- and 
double-junction systems^^'^^ as well as for tunnel junc- 
tion arrays^^'^". 

It is interesting to observe that both interaction correc- 
tion to the conductance and the shot noise spectrum^ are 
proportional to the same parameter /3 (0. This obser- 
vation illustrates a close relation between quantum noise 
and interaction effects in coherent conductorsSJ^. Pro- 
ceeding further along these lines one can investigate the 
effect of electron-electron interactions on current noise^-. 
The interaction correction to the Nyquist noise was again 
found to scale with the parameter (in accordance with 
the fluctuation-dissipation theorem) while the same cor- 
rection to the shot noise turned out to be proportional 
to the third cumulant of the current operator. It was 
conjecturedi^ that the same rule should apply for higher 
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cumulants as well, i.e. the interaction correction to the 
n-th cumulant should be proportional to the (n + l)-th 
cumulant for all values of n. A general pro of of th is con- 
jecture was very recently provided in Refs. \2Y22. 

Throughout the analysis^'^'^4 it was assumed that in 
the absence of interactions the conductor is described by 
an energy independent (though otherwise general) scat- 
tering matrix^"^. While in a number of important cases 
the above assumption indeed applies, in various other 
physical situations it turns out to be insufficient. There- 
fore it would be highly desirable to develop a generaliza- 
tion of the path integral techniqu o^i^° to conductors de- 
scribed by energy dependent scattering amplitudes. This 
generalization is the primary goal of our present paper. 

In physical terms our analysis should now effectively 
account for internal dynamics of coherent conductors or, 
in other words, for a finite dwell time of electrons. This 
effect should be combined with that of electron-electron 
interactions. Below we will accomplish this program for 
an important and widely studied class of conductors - 
the so-called quantum dots^'^'*. 

The structure of the paper and our main results are as 
follows. 

In Sec. 2 we describe our general real time path in- 
tegral formalism and derive the effective action of inter- 

m 



acting quantum dots, Eqs. H17I22|I of our paper. The ef- 
fect of electron-electron interactions is treated in a stan- 
dard manner by the Hubbard-Stratonovich decoupling 
of the Coulomb term in the Hamiltonian and reducing 
the problem to that of an electron interacting with the 
fluctuating quantum electromagnetic fleld defined on the 
Keldysh contour. In order to handle the fermionic part 
of the problem we combine our path integral analysis 
with the scattering matrix approach. For the model of a 
quantum dot adopted here the latter approach allows to 
exactly integrate out all the electron paths and express 
the effective action only in terms of the fluctuating fields 
which are then treated within an effective quasiclassical 
approximation suitable for highly conducting quantum 
dots with 5^1. 

In Sec. 3 we derive a general expression for the current 
noise in chaotic quantum dots as a function of frequency, 
voltage and temperature. Although this expression itself 
does not include interactions, it turns out to be very use- 
ful for better understanding of the relation between shot 
noise and interaction effects in mesoscopic conductors. 

Sec. 4 is devoted to a detailed description of the 
current-voltage characteristics of interacting quantum 
dots. In the most simple voltage-biased limit the cor- 
responding general expressions are presented in Sec. 4 A 
by Eqs. H53|l - H61() . Further analysis of these expressions 
for chaotic quantum dots in the leading non-trivial order 
in 1/g is carried out in Sec. 4B and 4C. A more general 
case of an arbitrary external impedance is considered in 
Sec. 4D. It is demonstrated that, as one goes away from 
the voltage-biased limit, the interaction correction to the 
current gets substantially modified and new regimes be- 
come possible. 



Further discussion of our results and their compari- 
son to several recent experiments can be found in Sec. 
5. In Appendices we present various technical details of 
our derivation of the effective action (Appendices A and 
B), details of our averaging procedure (Appendix C) and 
general expressions for the current in interacting quan- 
tum dots (Appendix D). 



II. EFFECTIVE ACTION AND CURRENT 
OPERATOR 

A. General formalism 

Let us consider a system of interacting electrons de- 
scribed by the Hamiltonian 



V2 



lldr 



dr' {r)^l (r') V^.' (r>, (r) . (2) 



Here the operator (V'cr (''')) creates (annihilates) an 

electron with the coordinate r and the spin projection a. 
In Eq. Q the summation over a is assumed. The time 
evolution of the density matrix of the whole system p{t) 
is determined by the equation 



p{t) = T exp [-iHt] p(0)t exp [iHt] 



(3) 



Applying the Hubbard-Stratonovich transformation we 
rewrite this equation in the form 



P{t) = 
where 



/ VVj Te-*^* ''*'-f^i(*')p(0)te*^' dt' H^(t')^^s^, 



/W.e'S 



. , „ Jo J 



(4) 



(5) 



is the electromagnetic contribution to the action, the 
symbols T and T imply ordering respectively in the for- 
ward and backward time directions and 



Hj = J dr^iir) 



2m 



- ^JL + U{r)-eV,{t\r) 



As usually, the expectation value of any physical op- 
erator is obtained by taking the trace of this opera- 
tor multiplied by p{t). Here we are mainly interested 
in the expectation value of the current operator. This 
operator, which we will specify later, can be obtained 
by taking the functional derivatives of the combination 



T(T)exp ±i£dt'Hi^2{t') 



with respect to the fluctu- 



ating flelds Vi,2{t). For this reason there is no need to 
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deal with the whole expression Q , the kernel of the cur- 
rent operator can be obtained from the effective action 
S, which is defined by the equation 



AS 



Tr ^Te~'fo ''*'-«"i(*')p(0)te*.''o dt'H2it')y 



(6) 



The effective action 5* can be expressed through the 
Keldysh Green function matrix G 



and the matrix Gy^ y^ satisfies the equation 

+ ^ + A' - Uir) + eV+it, r)^ 1 + SG'^ 

X Gv,M rr') - aj{t - t')5{r - r'), (8) 

Here we defined 5G^^ = {eV^ {t' ,r)/2)az (where cTz is 
one of the Pauli matrices) and introduced symmetric and 
anti-symmetric combinations of the fluctuating fields: 



(7) 



{Vi+V2)/2, V- =Vi-V2. 



(9) 



As we have already discussed, in the interesting for us 
limit of large dot conductances fluctuations of the fields 
Vi^2 remain relatively small. In this case one can expand 
the exact effective action in V^^ keeping only the first and 
the second orders. Then one finds 



iS — iSc 



2TT[Gy+SG-^] - Tr[(Gv+5G-i)2], (10) 
H) with 6G~^ = 0. 



where Gy+ is the solution of Eq 
From (IIOII we obtain 



;^ dt' J dPr{Gy+^^i{t' -Q,t',r,r) 



+ Gy+^22it' + 0,t',r,r))V-{t',r) 

-e^ f dt' f dt" f d^r'd^r"Gv+^i2{t',t",r',r") 
Jo Jo J 

X V-{t",r")Gy^Mt"^t'^r'\r'))V-{t\r'). (11) 

Note that at this stage it is important not to expand in 
keeping the exact nonlinear dependence of 110(1 on 
this field. 

The components of the Green function Gv+ can be ex- 
pressed through the initial single particle electron density 
matrix pQ and the single particle evolution operator . 
These are the matrices in the channel space denoted by 
a hat here and below. The operator Lf^ is evaluated in 
the Appendix A. We have 



Gy+,ll(il,t2) = ~ie{ti-t2W (ti,<2) 

+ iJ7'^^(ii,0)/5o?7^+(0,t2), 
Gv+M^iM) = -^0{t2-tl)U^^{tl,t2) 

+ iij^^{ti,{))pS^+{Q,h), 
Gy+42(ii,t2) = iiJ^\ti,Q)poiJ'^\Q,t2), 
Gy+,2i(ii,t2) = ~iU^* {ti,{))[i^ po]U^\Q,t2) 112) 



Vl= (pi_/e 



Vr= cpr/g 



Zs(a))/2 



Zs(a))/2 




FIG. 1: Schematics of the quantum dot. 



Here for future purposes we have defined 



^±(f ) = / dt"eV^{t"). 
Jo 



B. Effective action 



(13) 



The above analysis is rather general and can be ap- 
plied to a variety of mesoscopic structures with disorder 
and interactions. Our primary goal here is to consider 
electron transport through a system which contains an 
interacting quantum dot. This system is schematically 
displayed in Fig. 1. 

A quantum dot can be viewed as an island in-between 
two barriers connected to a voltage source via metallic 
leads with an arbitrary impedance Z{uj). Electrons can 
enter the dot through one of the barriers, spend some 
time there propagating between the barriers and being 
reflected, and finally leave the dot through another bar- 
rier. Details of the electron motion inside the dot will 
not be important for us: Electrons can either propagate 
ballistically from one barrier to the other or suffer addi- 
tional scattering inside the dot, e.g. from the outer walls 
or otherwise. 

Electron-electron interactions are taken into account 
by means of the effective capacitance model. Here we 
introduce capacitances of the left and right barriers Cl^r 
and the gate capacitance Cg . Making use of this model 
we can specify the expression for the effective action fur- 
ther. Namely, we will assume that strong changes of the 
fluctuating voltage fields in space are allowed only in 
the vicinity of the barriers where they suffer jumps 
and (see Fig. 1). Additional voltage drop inside the 
dot is neglected, i.e. there only the time dependence of 
the fields is taken into account. In the leads these 
fields are assumed to vary slowly in space. 
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Let us introduce the fluctuating potentials of the left 
and right leads VJ^ and of the dot Vj^ . Then we define 

the (time dependent) voltages across the barriers, = 
— and = — , and find the electromag- 
netic potentials in the rest of our system by minimizing 
the action . Afterwards we substitute the result back 
into Eq. and arrive at the expression which depends 

only on the phase jumps ^(i') — Jq dt" eV^ j^{t"). In 
this way two contributions to the total effective action 
can be identified, 

iS" = iS'dot + iS'ext- (14) 

Here the term iS'dot comes from the last two terms of 
Eq. where the space integrals run over the inner 



part of the dot and the barrier area. The second con- 
tribution zS'cxt comes from the integrals over the remote 
parts of the system (leads). Making use of a slow co- 
ordinate dependence of the fields in the leads one 
can expand the action iS'oxt m these fields and keep the 
first and the second order terms in this expansion. We 
arrive at the standard form of the action^^ describing 
fluctuations produced by an arbitrary environment with 
a linear impedance. Above we introduced the two leads 
with identical impedances Zs{u))/2. The fluctuating part 
of the voltage across the left lead is ((/jJ — )/e, the 
corresponding value for the right lead is ((/3^ -I- (p'^)/e. 
The corresponding contributions to the action need to 
be added and we arrive at the following expression 



'iScxt — 



6 Jo 



- 5 / dhdh Zs\h - t2) [{^lih) - ^;{h)){iptih) - ipj{t2)) + i^^ih) + ^;ih)){ip+{t2) + ^+(t2))] 
6 JO 



1 

e .10 



dhdt2 Zsih - t2) [{ipl{tl) - ip-{h)){^l{t2) - ip-{t2)) + (^«(tl) + ip~{h)){^-^{t2) + ip-{t2))\ 



2tt Zsiu^y 



27r 



2T 
V 



Re 



Zs{ijj) 



(15) 



The term iSam is evaluated within the capacitance model 
and included in the above expression for i^oxt H15|l . 

Now let us specify the effective action of the quantum 
dot. In what follows we will assume that the left (right) 
lead contains Nl {Nr) conducting channels. Every chan- 
nel is characterized by the velocity Vn (which may also be 
different in different leads). In the absence of electron- 
electron interactions quantum transport through the dot 
can be fully described by the scattering matrijii*^. As 
we have already pointed out in Sec. 1 a convenient way 
to include interaction effects is to combine the scatter- 
ing matrix approach with the path integral technique. 
As compared to Refs. "Pgl-f^ here we are dealing with 
a more complicated situation because the scattering ma- 
trix of the dot can now depend on the energy E of in- 
coming electrons. Therefore, the action^' cannot be 
directly used and a proper generalization is required for 
our problem. 

Let us define the scattering matrix describing electron 
transport through our quantum dot in the absence of 
interactions: 

For our derivation it is also convenient to perform the 



Fourier transformation of the S*— matrix. We introduce 



and the functions f(t), f'(t), f(t), t'{t) defined analo- 
gously. In accordance with the causality principle one 
has S{t < 0) = S\t < 0) = 0. 

It turns out that one can find an explicit expression 
for the evolution operator in terms of the scatter- 
ing matrix, construct the Green functions (|12|l and then 
derive the effective action of the quantum dot Sdot from 
Eq. (|ll|l . This program is carried out in Appendices A 
and B. The final result can be expressed in the form 



iS, 



dot 



iSn — St 



(17) 



As usually, the part Sr describes dissipative effects. It 
reads (see Appendix B): 



iSr — ~2i dz dxdy tr{[S{z — y)(p {z)d{z ~ x) 
Jo Jo 

-SHz-y)0-{z)S{z-xMy,x)}. (18) 
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where 



-lipliz) 
i^^(z) 



p{y,^)=P,{y-x)\ ^ (19) 



(1 is the unity matrix), and 



dE 6*-^"= 1 



27r 1 + e^^/^ 2 



2 sinh[7rTx] 



(20) 



is the equiUbrium density matrix of non-interacting elec- 
trons. 

The second term Sj in (|17|l accounts for quantum 
noise. This term is also evaluated in Appendix B. As 
a result we obtain 



Si = dxidx2 f dyidy2 f dzidz2 
Jq Jo Jq 

X tr{ [S{xi - zi)(p^{xi)S{xi ~ yi) 
-S^xi - zi)(p^{xi)S{xi -yi)]p{y2,yi) 

X [S{X2 - y2)^^ {x2)Six2 - Z2) 

-S\X2 - y2)V~ {X2)S{X2 - Z2)\h{zi,Z2)]i21) 

Here we defined 

Cgi[(/)+(2;i)-ip+(z2)] \ 

and ft.o(^) = 5{z) — pq(z). We also note that in the long 
time limit the integration over x, y in Eq. H18|l and over 
2^1,2, 2/1,2 in Eq. H21|l can be extended to the interval 

(—00,-1-00). 

The expressions (|17I22|) for the effective action of the 
quantum dot represent the main technical result of our 
paper. This effective action is defined by essentially 
nonlocal in time expressions which account for electron- 
electron interaction effects in the presence of a non-zero 
dwell time to of electrons in the quantum dot. Should 
T]j be much shorter than any other relevant time scale 
in our problem, the time dependence of the scattering 
matrices S{t) and S'^{t) can be approximated by the S- 
function, S{t) oc S~^(t) oc 6(t), in which case the action 
(|17l22f) reduces to one derived in Ref. (9t 



C. Averaging of the action 

If one is not interested in mesoscopic fluctuations of the 
effective action one can simplify the above expressions 
by averaging Eqs. H18|) and (|21|l over energy intervals 
exceeding the dot level spacing 6. Consider first the term 
Sr lfTH|) . Let us illustrate the main idea by treating the 
average tT{t'' {x)t{y)). We have 



tr(it(a;)%)) 



tr{i\Ei)i{E2))e^'''-'^^y.{2i) 



In a broad interval of energies the average ti:{t'^ {Ei)i{E2)) 
should depend only the energy difference Ei—E2- Making 
use of this observation let us define the function u^^ — 
{E)i{E + ui)) and its Fourier transform URL{t) — 
I l^u^^c^*'^*. The function u^^ satisfies the property 
""-f; = ' therefore URL{t) is real. Other averages are 
defined analogously. We find 

iv{[5{x)5{y)l - f\x)f{y)]) = 5{x ~ y)uLL{y), 

iT{i'\x)i\y)) = S{x ~ y)uLR{y), 

tr{[S{x)5{y)i~?\x)f'{y)]) = 5{x - y)uRR{y), 

tT{i\x)i{y)) = S{x^y)uRL{y).{24) 

Averaging of Eq. lfTS|l with the aid of if^ yields 



E 

i,j=L,R Q 



We note that - unlike the non-averaged action Sr 
(|18|l - is bilinear in both ip^ and (p~ . 

Now let us average the term Sj. We first notice that, 
since this term is already quadratic (p~ , in the averaged 
version of Si one can neglect fluctuations of the phases 
((f^^ and set f'^^nit) = eVL,Rt. After that the voltages 
Vl,r can be absorbed as energy shifts of the Fermi dis- 
tribution functions in the leads. Averaging of the com- 
ponents of the scattering matrix entering into Eq. (|21|l is 
carried out as above. We first proceed to the energy rep- 
resentation. Then we will get a sum of terms containing 
combinations of a similar structure, such as, e.g.. 



dE 
2^ 



f{E + LO~eVLKl- fiE-eVL)) 



xtr([i - fHE)f{E + uj)][i - f\E + uj)f{E)]) ,(26) 

where f{E) = (1 -t- exp(i?/T))~^ is the Fermi function. 
Since the averages should not depend on E, in all these 
combinations one can integrate over this variable. Col- 
lecting all terms we arrive at the final result 



ST = V / dxi f dx2 
, i Jq Jq 



Here the kernels v are defined as Vij (t) — J v'^ e 
where 



— coth — Rc 
27r 2T 



cj — eV , Lu — eV 
coth ■ 



An 



UJ 



eV UJ - 
■ coth ■ 



-eV 



2T 



UJ U) 

2^^°*^2T 



An 

V =^Vl + Vr and 

tT(f{E)r''{E)i'{E + uj)i'\E + w)). 



2T 



-tv{f{E)t^{E)f'{E + uj)t' {E + uj)), 
-tr(i{E)f'' {E)i'{E + uj)?\e + uj)) , 
iT{i{E)i\E)f' {E + uj)?\e + uj)). (29) 
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This concludes our derivation of the effective action for 
interacting quantum dots. 



D. Current operator 

In order to complete our general analysis let us define 
the kernel of the current operator for our problem. One 
can choose calculating the current either in the left or in 
the right junction, obviously in the stationary limit the 
result should remain the same in both cases. One can 
also use a symmetrized version of the current operator. 
One finds 



where 



1 / . SjtSn - Sj) . SjiSn - Sj) 

— I ~ZC " — IG 

2 V S^lit) 5^-j,{t) 



ie S{iSii-Si) 



2 



(31) 



Since here we are interested in a stationary situation, 
the result for the current should not depend on time. 
Therefore one can choose an arbitrary value of t at which 
the above functional derivative is taken, with the only 
requirement that t should be sufficiently large. Inserting 
Eqs. 1)18(1 and ((21|l into and defining the matrix 



A 



1 

-i 



we obtain 



lit, = Jo(i,^+)+(5/(i,^±), 



(32) 



where 



Io{t,f^) = e j dxdytv^[6{t- x)k5{t~y) 

-S\t-x)kS{t~ymx,y)], (33) 



and 



"V'r— const 



2" 



dx 



dyidy2 



dzidz2 Tr 



{ [S{t - zM5{t - y,) - S^{t - z,)AS{t - yi)]/5(y2, yi) 



X \i(x - y2)(p (x)8{x - Z2) - S^x - y2)(p {x)S{x - Z2)] h{zi,Z2)^ 



le 



dx dyidy2 dzidz2 Ti l[5{x ~ zi)(p {x)S{x - yi) - {x - zi)ip {x)S{x - yi)]p{y2,yi) 

^ -IQ JO JO ^ 



X [6{t - y2)k5{t - Z2) - S\t - y2)kS{t - 7.2)] h{zi, Z2)} . 



(34) 



r 



III. NOISE OF A QUANTUM DOT 

The general expression for the term Sj (121(1 in the ac- 
tion enables one to easily evaluate the current correlators 
at the barriers in the absence of interaction: 



S^,{io) = j dte''^\SI,{t)SIj{0) + SIjiO)dh{t)), (35) 
where i,j = L, R. One finds 



Sij{uj) = 2e^ 



dte" 



(36) 



In this way one recovers the well known general expres- 
sion for the noise in terms of the scattering amplitudes^,. 
Averaging the result over mesoscopic fluctuations we get 

2e2 



)w coth 



2T 



X (iu; - eV) coth ^^-^ + (c + eV) coth ^^^\s7) 



2T 



2T 



Following the standard procedure^^ let us express the 
scattering matrix of the dot in the form 



S{E) = R + f'[l - U{E)R']-^U{E)f, 
where the energy independent matrices 



and 



R = 



fL 

f' 



tL 
i' 



R' = i ''^ ^ 
' ffl 



4 
iji 



(38) 



(39) 



(40) 



account for scattering properties of the left and right 
barriers while U{E) is the unitary matrix which effec- 
tively describes scattering in the "internal" part of the 
dot. Further analysis requires specifying the model for 
our quantum dot. Here we will mainly address chaotic 
quantum dots in which case U{E) belongs to the circular 
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ensemble. The averages u'f and w,*-? are evaluated with 
the aid of the diagrammatic techniqueSSiSS. We obtain 



LL,BR 



,LR 



9L,R 

2 



gl.R 



2g{\ — ibJTr)) ' 

9L9R 

2g(l — iujtd) 



(41) 



The details of the derivation are presented in Appendix 
C. Here we defined the dimensionless conductances of 
the left and right barriers ql.r. = 2tr(i^ ^{E)tL^E.{E)) — 
2iT /c^Rlji, their sum g = + g^, and the electron 
dwell time in the quantum dot td = 'iir/gS. The general 
expressions for v^J can also be obtained with the aid of 
the results derived in Appendix C. One finds 



253 



9L 



'D 



1 + c^2r2 



9R 



+ 



9l9r {9l + 9Rf 



+ 9% 



9t9R 



2^4(1 +c.2r2) 

2 2 

9l9r 
2<?3(l+a;2r2) 



1 + wVI, 

Pr, 



+ 



+ 



gLgniii^TDgL + (1 + iujTB)gR] 
9l9r{{^ - ■i^TD)9L - ioJTDgR] 



Pi 



2g^{l+cu^rl) 
Here we introduced the Fane factors of the barriers 



(42) 



" tr(£,£t) ' ^ 



L,R. 



We note that Eq. H43|l agrees with the results for zero 
frequency noise of chaotic cavities previously derived in 
different limits by various author o^i^? . It is also satisfac- 
tory to observe that in the particular case of two diffusive 
conductors Pl — Pr = 1/3 Eq. again yields (3—1/3 
for any g^ and g^;. 

Our analysis also allows to generalize the results for 
the current noise in chaotic cavities to the case of finite 
frequencies. In order to find the noise correlator S{uj) it 
is in general necessary to account for the barrier capaci- 
tances. Assuming for simplicity Cg = 0, we obtain 



S(co) = 



where 



= \YLiio)\^SLLito) + \YRiLu)\^SRR{Lo) 
+ 2Re{YLiu)SLR{Lo)YRiio)), 



(45) 



>l(^) = 
YRiuj) = 



-iloCr+^4{u5^-u5^) 

-ioj{Cl + Cr) + f [u^L + uSR - u^R - uS^ 



r 



2 

-iljCl + — (m 



LL 



- u 



LR. 



MCl + Cr) + ^ (u^i + u«« - _ uS^) 



and u^J are defined in (|41|l . The above equations demon- 
strate that in a general case the noise spectrum (|45f) de- 
pends on frequency in a complicated manner. However, 
for fully symmetric quantum dots, i.e. for Cl = Cr, 
9l = 9r and Pl = Pr = /?, one find Yl{i^) = Yr{uj) = 
1/2. In this case all frequency dependent contributions 
contained in the functions wj' 1421) cancel out and Eq. 
(|45|l reduces to the standard form^ 



S{iu) = 2(l-/3)Goc^coth — 



-/3Go^(w±ey) coth^ 



•±eV 
2T 



(46) 



The functions and vl;^^ are recovered by interchang- 
ing the indices L ^ i? in the expressions (|42|l for v^-'^ 
and v^^. 

In the low frequency limit oj ^ one finds -D^^ = 



yRR 



^LR 



where S is the noise spectrum of the dot. 



^ ^ f 9l9r{9'^ - 9l9r) 



9l9rPl + 9l9rPr 



Rq \ 



.nth (slal , 9Lg%pL + gi9RPR\ 

R, '°*^2rlv 9' ^ 9' )■ ^ ' 

Here and below Rq = hje^ is the quantum resistance 
unit. 

In the leading approximation the dot conductance is 
Go = 9l9rJ gRq- Hence, for chaotic quantum dots the 
total Fano factor (3 (defined in a standard manner as a 
ratio between the shot noise spectrum and its Schottky 
value 2el) reads 



P- 



^ 9rPl+91Pb 

{gh + 9rY {gh + gRf 



gL9R 



(44) 



with the Fano factor /3 = (1 + /9)/4 in accordance with 
Eq. gH). 



and, hence, 5ll(0) = Srr{0) = 5, IV. 



CURRENT- VOLTAGE CHARACTERISTICS 
AND CONDUCTANCE 



Let us now turn to the interaction effects. In order to 
evaluate the current- voltage characteristics we should av- 
erage the above expressions for the phase-dependent cur- 
rent I{t,ip^) over the fluctuating phase fields 
ip^ . In addition to that - provided one is not interested 
in the effect of mesoscopic fluctuations ~ one can also 
average the result over such fluctuations. We proceed 
exactly in this order and first perform averaging over the 
phase fluctuations. 



A. Averaging over fluctuating phase fields 

Let us combine Eqs. ()32|l - (|34|) with (|30|l and carry 
out functional integration over ip^ . We notice that in 
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the interesting for us limit of large dot conductances it is 
parametrically justified to perform this integration with 
the averaged effective action 



nav 



instead of the exact one. This observation simplifies the 
whole procedure enormously because the action S^'^ is 
quadratic in ip^. Hence, the integrals become Gaussian 
and can be handled exactly. 

Since in all the integrals only linear combinations 
of the phases ip^ enter into the exponent, it is con- 
venient to first integrate out this variable. This in- 
tegration yields functional (5-functions. For example, 



one of the terms in Eq. H34(l contains the exponent 

ci[vtiy2)-vt(yi)+vU^2)-vU'^i)]^ in which case the averag- 
ing over </jJ fixes the variables (p^ in the form 

VjiT) ^ -KjL{zi - t) + KjL{z2 - t) - Kjuiyi - t) 
+ K,R{y2-T), j^L,R,g. (47) 

The Fourier transform of the functions Kij (t) reads 

^2 



^——A^iu;), 

— VjJ + 

where the matrix A{lo) has the form 



(48) 



A(lo) = 



^ —iujCh 



2 



Zs(^) 
2 



Zs(w) 



2 

Zs(w) 



(49) 



Other terms of Eq. H34() are treated analogously. 
The above expressions hold for an arbitrary external 
impedance Zs{lu). Further general expressions for the 
current are presented in Appendix D. 

Let us now focus our attention on the important limit 
of vanishing external impedance Zsito) — > 0. In this limit 
one finds 



K^^ = K 




where we defined 
= 



,LL 



- -iujCs + ^u{uj) 



= tT{i- S^{E)S{E + uj)). 



(50) 



(51) 



(52) 



and Cs = Cl + Cr + Cg . What remains is to integrate 
out the variable tp^ . Because of the obtained functional 
(5-functions this integration becomes trivial as well. One 
should just substitute the trajectories (|T7|l into Sf^ , into 
the imaginary part of 5ext and into the term SI. Tak- 
ing the unitarity of the S'— matrix into account, in the 



voltage-biased limit Zs{oj) — > one finds 

I{V)=Io{V) + 6IiV). (53) 

For the sake of convenience we present the expression for 
the term Io{V) in two equivalent forms: 



Io{V) = '^e j dxpo{x)T{x)e 



-F{x) ^^ieVlx _ ^-ieVRX'^ 



(27r)2 
xf7..e^"-(^) 



(54) 



Here we have defined T{x) = / |f tr[£t (£;){(£;)] e"'^^, 
and 



IfiE-eVL) \ 

^^=1 ifiE + eVn))' ^''^ 



where f{E) — 1/(1 -I- exp(i?/r)) is the Fermi function. 

The interaction correction to the current 5I{V) takes 
the form 



61 = -elm J dx J dyidy2 J dzidz2 J^{yi,y2, zi, Z2) {K{-yi) - K{-zi)) 



X tr{ [Sizi + x)k5{yi + x) - S\z^ + x)kS{yi + x)] Mvi - 2/2) [S{y2)S{z2) - ^^(^2)5(^2)] ho{z2 ~ zi)},(56) 

I 

where we have set In the limit Zs{lu) — > we have Re(iS'cxt) = and the 

T{yuy2,z,,Z2) = c-^iy^-y-)-P(^-^~^-2)'nvi-z.) 

X g--P(a2-Z2)+-F(yi-Z2) + -F(l/2-Zl) j-gy-J 
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function F{x) is simply equal to F{x) = evaluated 
for ip~[{T) = — <p^(t) = K{x — t) — K{—t). In this case 
we find 



F{x) 



V, 



duj 

LR 



-vS"-), (58) 



where the functions are defined in Eqs. H28|l - (|29|l . 

As we shall see below, for a wide range of parameters 
the function (|57|l in Eq. (|56l) can be approximated by 
unity ^ — > 1. Under this approximation one can con- 
veniently perform the Fourier transformation and reduce 
the expression H5t)|) to the form 



SI — elm 



dEduj 



K^tT{[l-S^iE)SiE + Lj)]hE+. 



X [S^ {E + uj)AS{E + uj)-S'' {E)AS{E)]pE } , (59) 

where He = 1 — pe- 

Employing the above general expressions one can also 
define the zero bias conductance of the dot. In doing so, 
we express the voltages Vl.b. in the form Vl,r = {V zt 
Vg)/2, where Vg is the effective potential of the dot tuned 
through the gate, and find the linear in V correction to 
the current. The voltage Vg can be removed by shifting 
the energy and performing a proper re-definition of the 
S"— matrix. Thus we find G = Gp + SG, where 

= -y/ ^Tr{i-S^E + co)AS{E + co)A} 



dl 

dE' 



(60) 



and 



SG = —Im 
2 



dEduj , , ,,df(E) 



(27r)2 '"^ " dE 

X tr{[S'+(£;)AS'(£;)A - S+{E + uj)AS{E + uj)A] 
•k[1- S+{E)S{E + uj)]]. (61) 

In Eq. (|61|l we again set = 1. This is sufficient in the 
perturbative regime to be studied below. 



B. Perturbation theory and renormalization of the 
S- matrix 

The above general expressions allow to conveniently 
proceed with the perturbation theory in the interac- 
tion. In order to obtain the first order correction to 
the conductance one should expand Gp (IHO)) in F{x) 
and combine linear in F(x) terms of this expansion 
with the interaction correction SG (|61|) . Then one finds 
G = Gl + SGi + SG2, where G ~ Gl is the Landauer 
conductance of the dot in the absence of interactions. 



dE 



df 



GL^~^^ir[l^S\E)AS[E)A]^, (62) 



and the two corrections 5Gi,2 read 
°2 r dEduj 



I , , E + u} 

K,,, coth tanh 

2T 2T 



df{E) 



iY{S\E + Lu)AS{E + iu)A - sHe)AS{E)A}^^, (63) 

Oil/ 



SG2 



Im 

2 



dEdw ^ E 

Ku, tanh 



d.f{E) 



{2n) 



2T dE 



X tY{S^E)AS{E)ASl{E)S{E + lu) 
-S''{E)AS{E + Lu)A}. (64) 

We note that in Eq. H63() the terms containing 
coth(a;/2r) come from the expansion of H6U|I while all 
the tanh-terms originate from Eq. (|61|) . 

It is easy to see that SGi = for any S{E) and at 
any temperature. Indeed, making a shift E + lo ^ E 
in the term involving the product S^E + u;)AS{E + uj)A 
one reduces the expression under the integral to the form 
ImKi^ B{E, Lu), where 



BiE,Lo) 



coth tanh — ) 

2T 2T J dE 

E + uj\ df{E) 



coth tanh 

2T 2T 



dE 



(65) 



Since f{E) is the Fermi function, one has B{E,uj) = 
B{E,—u!). At the same time, as follows from Eq. H51|l . 
ImK^ is an odd function of lo. Hence, the integral over lo 
vanishes identically, the term SGi drops out and only the 
combination H64|) needs to be taken into consideration. 

Let us make a shift S{E) SiE) + SS{E) in the 
Landauer formula 1)62(1 . Then the corresponding linear 
in SS{E) correction to the conductance becomes 



^Gl = ^tr[SSHE)AS{E)A 



- SHE)ASSiE)A]^. 



(66) 



Comparing this formula to the first order correction H64() 
we can choose SS{E) and SS^'lE) in the form: 



te.nh^±^[S{E + u;)~S{E)] 



tanh SiE)[S^{E - uj) - S^{E)]S{E) 



SS^ 



2T 

-f-K 

2i J 2tt 

E~uj 



tanh {S\E^u)^ {E)\ 



2T 



tanh ■ 



2T 



S\E)[S{E~u)- S{E)]S\E) 



(67) 



Although this choice is obviously not a unique one, we 
will stick to it for the reasons which will become clear 
below. Making use of the property K^^ = K*, it is 
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easy to check that the conditions {5S^{E)y = 6S{E) 
and SS''{E)S{E) + SHE)SS{E) = are identically ful- 
filled, i.e. the first order renormalization (|67|l preserves 
the unitarity of the S*— matrix. 

In order find the function K^j (151(1 for chaotic quantum 
dots one needs to perform the average of the product of 
two S*— matrices at different energies Neglect mg 

small localization corrections, and applying the formula 
31 we obtain 



ti{i- S^{E)S{E + uj)) 



'IbJTog 



2(1 — ILUTjj) 



Combining with (|51I52I) . we get 



(68) 



(69) 



and after the the Fourier transformation we arrive at the 
expression 



K(t) 



2tt 
9 



0{t) 



'D\ 



-t/r 



t 



{td + Tof 



To + TD 



(70) 



Here we have defined the classical i?C— time of the central 
island tq — 2-KC^/{e'^g) and r = tdToI{td + tq). For 
large quantum dots one typically has tq ^ td and, hence, 

T « Tq. 

Eq. H70|) allows to specify the necessary condition of 
applicability for our analysis. Let us recall that during 
our derivation we have expanded the effective action in 
powers of the fluctuating phase field tf~ . The least action 
condition fixes this field to be equal to ip~{t) = K{t). 
Since all the available time integrals are effectively cut at 
times exceeding TD^ it is sufficient to require lp~{td) — 
K{td) ^ 27r. Then from Eq. ifTUIl we conclude that our 
expansion of the exact effective action in powers of the 
field ip~{t) is justified in the "metallic" limit 



5> 1- 



(71) 



We also note that this condition is necessary but in gen- 
eral not a sufhcicnt one in order to justify the perturba- 
tive expansion in the interaction employed in this subsec- 
tion. An additional condition will be established below 
in Sec. 4C. 

Consider the most relevant physical limit td 3> tq. In 
this case one has w — ( —. , / ^ ■ ] ■ Since the 

integral in Eqs. (|67|l is taken over a wide range of en- 
ergies, the main logarithmic correction to the scattering 
matrix can be obtained if one replaces S{E -\- uj) by its 
energy independent average {S{E)) = R. Then one finds 

5S[E) = '^^^^1^[{S) - S{E) 
9 

-S{E){{S^)-S\E))S{E)\. (72) 

Eq. H72|) can be used to derive renormalization group 
(RG) equations for the energy dependent scattering ma- 
trix of a quantum dot. Following the standard procedure 



let us fix an infinitesimal energy interval between E and 
E — dE. From Eq. 1)72(1 one easily finds the correspond- 
ing correction dS{E) to the scattering matrix. Lowering 
the energy and repeating this procedure many times one 
arrives at the following RG equations 



dS{E) _ (S) - S{E){S^)S{E) 



d\n{l/ETQ) 



(73) 



where the effective conductance g is expressed via the 
rcnormalized scattering matrix at a given energy by 
means of the Landauer formula. Making use of Eqs. 1(381 - 
140(1 we define 
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= trf . 



(74) 



Hence, g itself gets rcnormalized and becomes energy de- 
pendent. Finally, substituting S{E) in the form 1(38(1 we 
observe that U{E) remains unchanged in the course of 
renormalization. Thus Eq. ((73() describes the renor- 
malization of the barrier transmissions connecting the 
chaotic cavity to the ideal leads. After simple manipula- 
tions and making use of Eq. 1(741) we arrive at the scaling 
equations for the transmission matrices T+T: 



dlii{ETo) 



trTtT 



(75) 



Eq. ((75() coincides with that recently derived in Ref. |23 
by means of a different approach. 

Let us emphasize again that the whole analysis of the 
present subsection is valid within the first order pertur- 
bation theory in the interaction. Hence, the above RG 
equation can be applied only as long as the rcnormalized 
conductance g{E) remains large. 



C. I-V curve in the voltage-biased limit 

Provided temperature and/or voltage exceed the level 
spacing in the dot and provided one is not interested in 
resolving subtle details of the I—Vcmve related to meso- 
scopic fluctuations, it is convenient to perform averaging 
of the above general expressions for the current over such 
fluctuations. The average values of the products of re- 
flection and transmission matrices depend only on us, not 
on E. Hence, one can integrate over E and derive the 
current-voltage characteristics of a quantum dot in the 
presence of interactions. In particular, the term Iq in 
Eq. ((54(1 acquires the standard Landauer form 



lo = -tY{iHE)i{E))V = GoV. 

TT 

In the case of chaotic quantum dots one finds 



(76) 



tr(i\E)i{E)) = 



9L9R^ 

2.9 



,9l9r{'^ - M + 9l9B.il - Pr) 



(77) 
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where v = —1, or 1/2 respectively for circular orthogo- 
nal, unitary and symplectic ensembles. The second term 
in Eq. iffTjl represents the weak localization (WL) correc- 
tion. For the model under consideration this correction 
is parametrically (~ l/ff) small as compared to the first 
term in (|77|l . Below we will concentrate on the interac- 
tion correction to the current 51, which is of the same 
order in 1/g, but nevertheless strongly exceeds the WL 
correction at sufficiently low energies. 

Let us consider the perturbative in the interaction 
regime. After averaging over mesoscopic fluctuations the 
interaction correction (|59|l takes the form 



51 



^ lu,{K^D{u^)) [u - eV) coth 



(78) 



where we defined 

D{uj) = iv([f+ {E + uj)i' {E + uj) - f+ {E)i' {E)] 

X [i'"" {E)f{E + Lj) + ?^ {E)i{E + uj)\){Td) 

All other contributions to the interaction correction van- 
ish upon averaging. An explicit evaluation of the function 
D{ijj) (|79|l is easily performed with the aid of Eq. ljC9|l . 
Here we only quote the result: 



D{u) 



gLgRigB.PL + gLPn) 



'2{gL + gnY (1 - ioJTD? 
Combining Eqs. l(7S|) . and (|Snj), we obtain 

eB f duj 1 / 1 1 
~2~ J 2^^\l+uj^T^ 
uj + eV 



(80) 



51 



1 + lu^tI 



[U! ■ 



eV) coth ■ 



{uj — eV) coth ■ 



eV' 



2T ' 2T 

where we defined 

B = gLgnigLPR + gRPL)/{gL + gRf- 
One can also rewrite Eq. (|81ll in the form 

TT 



,(81) 



(82) 



51 



dt- 



sinh TiTt 



1 



e"*/^°)sineT^t. 

(83) 



The integral can be evaluated analytically. We find 



51 = Im 

TT 



- + ieV 



* 1 



1 



27rTr 



— + ieV 

TD 



^ 1 



1 



eV \ 
2t:T ) 
. eV \ 



+ 1 



2t:Ttd 2'kT 



,(84) 



where 4'(x) is the digamma function. This is a complete 
expression for the first order interaction correction to the 
current in chaotic quantum dots with large conductances 
in the voltage-biased regime. 

We will now analyze the above general expression in 
various specific limits. At T — > the I—V curve reduces 
to the following simple form 



d]_ 
dV 



Go 



B 
R„ 



■In 



T^{l + {eVTDy)- 



(85) 



Provided the voltage is large eVr ^ 1 the interaction 
correction turns out to be small. In the limit tq <^ tjj 
one finds 



dl 
dV 



— Go — 



Bg'El 



(86) 



where Ec — /2Cy: is the charging energy. In the in- 
termediate range of voltages 1/td ^ eV <^ I/tq the 
interaction correction becomes logarithmic in V 



dl_ 
dV 



Go 



2B 



In(eyTo) 



(87) 



and at lower voltages and temperatures eV, T <C 1/td 
we obtain 



dl 
dV 



= Go- 



2B 
To 



In 



{eVf 



(88) 

The latter expression demonstrates that in the regime 
under consideration the linear conductance of the quan- 
tum dot remains non-zero down to zero temperature. 

Taking the limit V" ^ in the general results (|83I84() 
one can derive the linear conductance of the dot at arbi- 
trary temperatures. Evaluating the integrals, we obtain 

G = Go - ^ [L{Tt) - L{Ttd)] , (89) 

JXq 



where we defined 



L{x) = * 1 



2tix 



1 

2iix 



(90) 



and 7 ~ 0.577 is the Euler constant. As before, this 
general expression can be simplified further in various 
limits. At high temperatures T ^ gEc we find 



G = Go 



Bg j Ec 3C(3)r7 ( Ec 



Rq I 3T 2tt^ 



T 



(91) 



where C(3) ~ 1.202. In the regime 1/td < T < I/tq we 
again arrive at the logarithmic interaction correction to 
the conductance 



2B 

G - Go + — In (Tro) 
Ho 



(92) 



which crosses over to Eq. (|88|) as T becomes smaller than 
the inverse dwell time 1/td- 

The above results completely describe perturbative in- 
teraction correction to the current-voltage characteris- 
tics of highly conducting chaotic quantum dots in the 
voltage-biased regime. We also note that in certain limits 
these results reduce to ones obtained previousl y^i^^i^'^i'^° 
by means of different techniques. For instance, the log- 
arithmic behavior of the interaction correction (|92|) can 
easily be recovered from the RG approach— which we 
also discussed in Sec. 4B. On the other hand, some other 
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results, e.g., Eqs. H86I91() . cannot be obtained within this 
RG scheme. 

Under which conditions do the above perturbative re- 
sults remain valid? At high temperatures the validity of 
the result (|91|l is guaranteed by the inequality T ^ gEc- 
In this limit no additional requirement (|71|l is needed, 
i.e. Eq. H91|) applies also for g ^ 1. The same is true 
for the result H86|) derived in the limit of large voltages. 
On the other hand, at lower temperatures and voltages 
T, eV <^ gEc the condition (|7T|l should be satisfied en- 
suring that fluctuations of the phase cp^ (t) remain weak. 

An additional applicability condition for the above re- 
sults is obtained if one requires fluctuations of the phase 
to be small. Let us recall that in the course of 
our derivation we have approximated the function H57|l 
by unity. This approximation is appropriate as soon as 
the function F{t) — {{ip^{t) — (^+(0))^)/2 remains much 
smaller than one and fluctuations of the phase (p+ (t) can 
be disregarded. Hence, the condition 



Firp) ~ 1 



(93) 



defines yet one more time scale in our problem, Tp , which 
restricts the applicability of the above perturbative re- 
sults. 

In order to determine the time scale Tp it suffices to 
evaluate the function F{t) in the limit of low voltages 
V ^ 0. In this limit Eq. reduces to a simple expres- 
sion 

^Im(if^)coth^(l-costjt). (94) 
Combining this expression with Eq. H69|l one finds 



g[To + td) 



2k , uj coth^ 1-cosLjt 



1 +U^T 



2_2 



(95) 

where k = T^jirp) -|- tq)'^. In the most interesting limit 
T ^ the integral in Eq. (|95|) yields within the loga- 
rithmic accuracy F{t) ~ (2^/^) ln(t/T). Then from Eq. 
we obtain 



Tp ~ rexp(5/2K). 



(96) 



Thus fluctuations of the phase <y5+(t) can be neglected 
at all temperatures and voltages only provided the dwell 
time Tp) is much smaller than the parameter Tp (|95|l . 
More generally, the results presented in this subsection 
are correct under the condition 



max(T,ey,r^^) > 



(97) 



It is easy to check that in this case the interaction cor- 
rection in Eqs. H87|) . (|88|l and (|92|l remains small as com- 
pared to Gq. If the condition (|97(l is violated the system 
enters an essentially non-perturbative regime in which 
case more accurate analysis becomes necessary. This 
analysis is beyond the scope of the present paper. Let us 
also point out that, while the non-perturbative regime is 
important from a theoretical point of view, it does not 
seem to be of much practical relevance for quantum dots 
with g ^ 1 considered here. 



D. Effect of external impedance 

So far all our flnal results have been formulated for 
the case of ideally conducting external leads ^s(w) 0. 
In many experimental situations, however, this voltage- 
biased model is not appropriate since the impedance of 
external leads remains non-zero in the relevant frequency 
range. For this reason it is important to find out how the 
above results for the 1~V curve are modified in the case 
Zs{^) + ^. 

In order to answer this question we will make use of 
the general expressions for the current derived in Ap- 
pendix D. Averaging these expressions over mesoscopic 
fluctuations we arrive at the result 



J(F) =/o(F)+(5/tot(V^), 



(98) 



where /o(^) remains the same and is given by Eq. H7t)|l . 
while 5/tot(^) represents the total interaction correction 
which takes the form: 



2T 



(^ - eV) coth Im (i^)D,, (^), (99) 



where 



Drl{u) = lx{i'{E + uj)i\E)f(E^i^)r'^{E^u)), 

Dll{lu) = tr{r\E)i'{E)i'\E)f{E + Lu)), 

Drr{oj) = tr{i{E + Lj)i^E + uj)f'{E + uj)?\E)), 

Dlr{oj) - iY{i{E + oj)i\E)f'{E)r'\E)). (100) 

The above expressions are still rather cumbersome. For 
the sake of simplicity below we will analyze the case of 
symmetric quantum dots with gp ^ qr — g/2, (ip = 
13r = Rp = Rr = R and Cp = Cr = C. We will also 
assume the external impedance to be Ohmic ZgiiJ) = Rg 
at all relevant frequencies. In this case by virtue of Eq. 
U wc find 



Dpp = Drr = 
DpR = Drl 



gP , .9(1-2/3) 



g(3 



16 m{l~iujTD) mil-iuJTDY' 
.9(1 + 2/3) 5/3 



32(1 -iwTc) 32(l-iwTD)2 



101) 



From Eqs. gSJEH) we also obtain Kpp = Krr ^ Ks + 

K, KpR = Krp ^ Ks~ K, where 



K 



-iiOCn 



2 —iUTD 



{-iLu + 0){~2iu;C+^ + i 



iujRsCg /A R l — iujTo 

(102) 



To simplify the analysis further let us replace 1 — 
iwRsCg/A by 1 in the denominator of the expression for 
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K (|102() . This is appropriate either for Cg ^ or in the 
hmit Rs < 2R. Then we get 

<5/tot(^^) = 61 {V) + SIsiV), (103) 

where SI{V) is the interaction correction evaluated for 
Zs ^ and 6Is{V) represents an additional contribution 
due to the external shunt. The term SI{V) is defined by 
Eqs. ^81I83I84|I . whereas for SIsiV) one finds 



X 



(w + eV) coth 



oj + eV 



{u! — eV) coth 



lOJTj) 



(104) 



2T ' ' 2T 

This formula can also be transformed to the following 
expression 

/ dt- 



SIsiV) -- 
(3{1 



X 



4tt g + Ags 



sinh-" nTt 



■ sin eVt 



Tse 



-t/rs 



, (105) 



where we defined ts = RRsC / {Rs + 2R) . As before, this 
integral can be expressed in terms of the ^'-functions. 
For the sake of brevity we will omit the corresponding 
expressions here. 

Let us concentrate on the linear in voltage regime. 
Making use of the above general results one can easily 
determine the total linear conductance G of the dot in 
the presence of an external Ohmic shunt. It reads 



G — Go 



1 

td{L{Ttd) - L{Tts)) 



(3{L(Tt) - L{Ttd)) 



9 



Td — TS 



+ (/? + 1)L{Tts) 



(106) 



where the function L{x) was defined in Eq. H90|) . 

For large quantum dots one typically has td 3> tq > 
Ts- In this case at sufficiently low temperatures Ttd <ti 1 
the whole quantum dot can be considered as a single 
scatterer, and we find 



G — Gn 



2/? 



Rqg + 4g= 



7 + 1 + In 



1 



2ttTtd 



where 



Go — Go 



JL 

2Ra 



In 



TD 



,(107) 



(108) 



and [3 ~ (/3 + l)/4 is the effective total Fano factor of a 
symmetric chaotic quantum dot, see also Eq. 1441) . Eqs. 
p07|l . (|l()8|l smoothly match with the result^ derived for 
a single coherent scatterer with td 



< 



To- 



In the opposite high temperature limit Ttd S> 1 , albeit 
Ttq <C 1, we obtain 



G ~ Gq 



f3 



2Rq 9 + 4gs 
P 9 + 2.95 



■In^ 



Rq9 + 



Ts 

7 + 1 + In 



1 



27rTro 



(109) 



100 - 
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FIG. 2; Conductance of a symmetric quantum dot, as given 
by Eq. 11061 . The Fano factor of each barrier is chosen to 
be /? = 1/3, its conductance is 200/-Rq (hence Go = 100/7?^) 
and td/tq — 10*. Different curves correspond to different 
values of the shunt resistor: (a) Rs = 2.5 x 10~^ Rq, (b) 
Rs = 0.0025 i?g, (c) Rs = 0.005 i?,, (d) Rs = 0.025/?,, and 
(e) Rs — 0.25Rq. A crossover between the two logarithmic 
regimes is clearly observed at Ttd ~ 1- 



Let us compare the results H107(l and H109|l . In both 
temperature regimes we observe that the interaction cor- 
rection depends logarithmically on temperature, though 
with different prefactors in front of the logarithm. At 
high temperatures Ttd ^ 1 this prefactor is determined 
by a sum of two different contributions originating from 
the terms 61 and Sis- In this regime the two barriers 
of the quantum dot behave as independent ones. In the 
other regime Ttd < 1 the first logarithmic term sat- 
urates and does not depend on T anymore providing 
effective renormalization of the non-interacting conduc- 
tance (|108|) . On the contrary, the logarithmic tempera- 
ture dependence (|l()7|l in the second term survives down 
to exponentially small temperatures. In this regime the 
quantum dot behaves as a single coherent scatterer. Its 
internal structure becomes insignificant in this case and 
the interaction correction scales with the total Fano fac- 
tor /3. A crossover between the two logarithmic regimes 
is clearly observed in Fig. 2. It is also important to em- 
phasize that, while the logarithmic in T correction H107|l 
vanishes for Rs — > 0, in the opposite limit Rs 3> i? it 
becomes practically independent of the shunt resistance. 

The presence of two different logarithmic regimes can 
also be understood bearing in mind a close relation be- 
tween shot noise and interaction effects in mesoscopic 
conductors. At high temperatures Ttd 2> 1, i.e. in the 
regime of independent barriers, the noise at each of them 
determines the corresponding contribution to the inter- 
action correction. Hence, the latter is proportional to 
the parameter p. In the opposite limit Ttd <SC 1 the bar- 
riers are not anymore independent. In this regime the 
temperature dependent part of the interaction correction 
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should be related to the total shot noise of the quantum 
dot which, as we discussed in Sec. 3, depends on the 
parameter 13. The high frequency noise wr^i > 1 now 
only provides the renormalization of Gq which is again 
proportional to (3. 

We will continue the discussion of the above results 
and their experimental relevance in the next section. 



V. DISCUSSION AND CONCLUSIONS 

The analysis presented here provides a general theo- 
retical description of electron transport through disor- 
dered interacting quantum dots in the "metallic" limit of 
large conductances. Our formalism combines real time 
path-integral-based influence functional technique with 
the scattering matrix approach. The main technical re- 
sult of the present work is the derivation of the effec- 
tive action for quantum dots described by arbitrary en- 
ergy dependent scattering matrices. Due to this partic- 
ular feature our technique allows to investigate quantum 
transport of interacting electrons in a very wide class of 
conductors which includes, e.g., structures with resonant 
transmission and many others. This formalism is a di- 
rect generalization of our earlier approachSiili which em- 
braced structures described by energy independent scat- 
tering matrices and did not account for internal dynamics 
of disordered conductors. 

Although our analysis is mainly aimed at interaction 
effects, the first and immediate result of our derivation 
is the expression for frequency-dependent current noise 
in non-interacting quantum dots. In fully symmetric 
dots the shot noise spectrum turns out to be defined by 
the standard expression (|46|l with the total Fano factor 
(3 = (/3 -I- l)/4. For asymmetric quantum dots the corre- 
sponding expression becomes far more complicated. 

Turning to interaction effects, let us briefly summarize 
our main results for the linear conductance G of chaotic 
quantum dots (characterized by dimensionless conduc- 
tance 5^1 and charging energy Ec) connected to the 
voltage source via an Ohmic resistor Rs- The conduc- 
tance G is expressed in a general form 

G = Go + (5G(T), 

where Go is the conductance of a non-interacting dot and 
5G is the interaction correction. This correction is neg- 
ative, i.e. electron-electron interactions tend to suppress 
the conductance of quantum dots. 

For the situation under consideration one typically has 
TD » To ~ l/gEc > Ts = RsC/2{l + GoRs)- For sim- 
plicity we quote the results for symmetric quantum dots. 
At high temperatures Ttq 3> 1 in the leading approxi- 
mation one finds 



SG/Go 



3T ' 



(110) 



where f3 is the Fano factor of a single barrier and 
the parameter x is defined below. At lower temperatures 



Ttq <^ 1 the power law dependence (|110|l crosses over 
into the logarithmic one 



SG^^^HTro). 



(Ill) 



The parameter x iii Eqs. Hll()|) . (|lll|l depends on the 
relation between Rs and the dot resistance Rd = 1/Go. 
For Rs — > one has x = Ij while in the opposite limit 
Rs 2> Rd this parameter is x = 2. Thus, in both limits 
pi()|l and Ulllfl the magnitudes of the interaction cor- 
rection SG evaluated in the current- and voltage-biased 
regimes differ by the factor 2. The magnitude of SG is 
smaller in the latter regime because voltage fluctuations 
across the dot are suppressed in this case. 

Eq. Hlll|) applies down to temperatures T ~ I/tjj. 
At lower temperatures another logarithmic regime sets 
in and the difference between current- and voltage-biased 
situations becomes more dramatic. In the limit Ttd <C 1 
one finds 



G — Go 



2(3 



1 



Rq 1 



4.Rd/Rs 



MTtd), (112) 



where 



Go ~ Go - TTTT In — . 



2Ra 



To 



(113) 



Eq. (|113|l describes an effective renormalization of Go 
by electron-electron interactions. This renormalization 
is again twice as big in the current-biased regime as it is 
in the voltage-biased one. More importantly, as follows 
from Eq. I|112l) , for any non-zero Rs no saturation of the 
the logarithmic dependence of the interaction correction 
on temperature is expected at Tt^ < 1. We also observe 
that at Ttjj ^ 1 the interaction correction SG always 
scales with P and, hence, vanishes completely for fully 
transparent barriers /3 — > 0. By contrast, at lower tem- 
peratures Ttd ^ 1 only the T-independent renormaliza- 
tion of Go disappears in this limit, while the correction 
(|112|l scales with f3 and remains finite even for f3 — 0. 
This fact illustrates a direct relation between shot noise 
and interaction effects in mesoscopic conductors, cf. Eqs. 

and (|112|l . We also observe that in the hmit /3 — > 
the conductance H112|l becomes completely independent 
of the charging energy Ec- 

Similar results can be obtained for a non-linear / — 
V characteristics of our system. In particular, at T — > 
the corresponding expressions are reproduced if one 
substitutes G dl/dV and T ^ eT^ in Eqs. |)110I113|) . 

Both logarithmic regimes in the dependence of the 
interaction correction on temperature and voltage dis- 
cussed above were observed in many experiments on var- 
ious mesoscopic structures. Here we will briefly discuss 
only a few examples. In the experiment''^ this depen- 
dence was found in short diffusive metallic bridges at 
temperatures (voltages) T < IK. No saturation was ob- 
served down to the lowest temperature < 100 mK. Es- 
timating the Thouless energy 1/t^) for the parame- 
ters of Ref. i31i one arrives at the value of order few 
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Kelvin. Hence, this experiment was performed in the 
regime Ttd < 1 (|TT^ . 

Recently it was argued^ that in order to explain 
the experimental results-^ it is necessary to assume 
the presence of insulating barriers at interfaces between 
the bridge and the reservoirs. We note, however, that 
this assumption can hardly be justified because for the 
samples^^ it would inevitably lead to values of g 10-^ 100 
times smaller than actually measured. The presence 
of insulating barriers can also be ruled out on the ex- 
perimental grounds22i. This controversy is eliminated if 
one recalls that the calculation's was performed within 
the voltage-biased model Rs — *■ 0. For this reason 
the authors^ii have overlooked the logarithmic depen- 
dence (I112|l and attributed the observations^^ to the 
regime Ttd ^ 1, hence, requiring much smaller values of 
1/tu < 100 mK. Since the voltage-biased model does not 
strictly apply to the experiments^^, no such requirement 
is actually needed. 

Note, that in some other experiments, e.g., ones with 
multiwalled carbon nanotubesSiSliS^, the voltage-biased 
model appears to be applicable. Both the logarith- 
mic temperature dependence^ of SG at larger T and 
its saturation at temperatures of order few Kelvin were 
clearly observed. The latter temperature range is con- 
sistent with the estimates of the parameter 1/tjj for the 
nanotubeo'^^i^'^1^^. 

More recently, the logarithmic temperature depen- 
dence of the conductance was observed in strongly disor- 
dered multiwalled carbon nanotubes^^''^®. No saturation 
of this dependence was found down to the lowest mea- 
surement temperature T ~ 100 mK. Although estimates 
of the parameter 1/td for the samples^S*^ yield values 
below 100 mK, it is not easy to interpret these observa- 
tions in terms of the high temperature logarithmic regime 
(|lll|l because Lt ~ y^D/T {D is the diffusion coeffi- 
cient) remains shorter than the length of the nanotube 
L at all relevant temperatures T > 0.1 K. In this situa- 
tion our zero-dimensional description cannot be applied 
to the nanotube as a whole. 

An alternative explanation is based on viewing the 
nanotube as a chain of N connected in series statisti- 
cally independent segments of the length Ls « L/N. For 
Lt ^ Ls each of the N segments should behave as a zero- 
dimensional coherent scatterer shunted by an external 
impedance effectively produced by the remaining A'^ — 1 
segments. This scenario applies provided the scale Ls 
exceeds an effective transversal dimension of the system 
Ltr- In the opposite case < Ltr a slightly more gen- 
eral picture of an array of M x A^ independent scatterers 
can be introduced, where the number Ad depends on the 
ratio Ltr/Ls- Applying the result \iVl\ for each of the 
segments and assuming A^ ^ 1, at temperatures Ttd ^ 1 
one finds the conductance of the system in the form^S. 

IB M 

G{T)^Go + f-HTrn\ (114) 
where Gq represents its conductance at T ^ I/td, M = 1 



for Ls > ^tr and M > 1 otherwise. For diffusive nan- 
otubes one has /3 « 1/3 and td ^ L^/D. The crucial 
point of this scenario is that at sufficiently low T the 
scale Ls does not depend on temperature and is set by 
interactions. 

Note that a similar effect was also observed in disor- 
dered metallic wires^S. Clear deviations from the stan- 
dard Id behavior of the interaction correction 5G oc 
1/Vt were reported at temperatures T ^ I/ti^, where 
is the low temperature dephasing time measured in 
the same experiment. For the parameters^ one can ver- 
ify that both the observed magnitude and temperature 
dependence of the interaction correction are consistent 
with Eq. H114() . where Ls ^ ^ \/L>t^ and M = 1. 

Turning again to the data^S, let us recall that in this 
case the observed values of were considerably smaller 
than the nanotube circumference ird. Therefore we can 
assume Lg < ird and define^ M ~ Tid/Ls- Then we 
obtain M/N ~ nd/L, i.e. Ls drops out and for (3—1/3 
the prefactor in front of the logarithmic term in Eq. H114|) 
becomes ~ 2.1d/RqL. This universal value fits well with 
the observations^ at T < 10 K. Finally, we note that 
the latter inequality is also consistent with the condition 
Tt^p ^ 1 because the measured values of were found to 

"^^ sek. Hence, the assumption 



^ 10- 

T^p appears to be supported by the 



be in the range r^p 
Ls ^ L^ and td ^ 
data^^ as well. 

Further experimental and theoretical investigations of 
the effect of interactions on low temperature electron 
transport in disordered conductors are warranted. 
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of their data and experiments2ii2Li£. This work is part 
of the Kompetenznetz "Funktionelle Nanostructuren" 
supported by the Landestiftung Baden- Wiirttemberg 
gGmbH. 



APPENDIX A: EVOLUTION OPERATOR 

For non-interacting systems the relation between the 
S"— matrix and the Green functions, or, which is the same, 
the evolution operator is known since the work by Fisher 
and Lee^. Here we use the inverse form of this rela- 
tion, namely we express the evolution operator in terms 
of the S*— matrix. We introduce the components of the 
evolution operator Uij{t, 0; r, r'), where i,j label the left 
(Z) and the right (r) leads as well as the dot (d). The 
component Uu reads 



1 r dE 



-iEt+iE(r/v'^-r' /v'l^) 



X [(Sa/3 + e{r)e{-r'){r„p{E) - S^fj) 
+ ei-r)9{r'){r+yE)-5ap)]. 
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Here a, (3 are the channel indices, is the velocity in the 
a— th channel of j— th lead. Three other components of 
the evolution operator Urr, Uir and Uri are expressed 
analogously. It is convenient to introduce the "flight 
times" in the a— th channel x — r/v^^, y = r'/v^ and 
express the wave functions and the evolution operator 
via these new variables. The velocities then disap- 
pear, and we can conveniently rewrite the above relation 
in the matrix form 

U{t,0;x,y)^iS{t-x + y) 

+ 0{x)e{-y)[S{t -x + y)- id{t ~x + y)] 

+ 0i-x)e{y)[S+{-t + x-y)- idit ~x + y)]. (Al) 

Note that here we present U in the form of a 2 x 2 (rather 
than 3x3) matrix leaving out the components describing 
transitions (e.g. from the dot to the dot and some others) 
which turn out to be unimportant for our analysis. We 
will return to this point in Appendix B. 

The operator 1A1|I has the following properties: (i) 
?7t(t,0) = t/(0,i), and (ii) lim U^t, 0)U{t, 0) ^ 1. At 

t — >+oo 

any finite time t this operator does not obey the unitarity 
condition. This is quite natural because the components 
describing the evolution in the dot are missing in Eq. 
(jAlj) . However the operator U{t, 0; x, y) becomes unitary 
in the limit f — > oo. 

Let us now introduce the fluctuating potential V{r,t) 
and, as we have already discussed in Sec. 3, assume that 
this potential is spatially constant both inside the metal- 
lic leads and inside the dot, and it suffers jumps at the 
junctions between the dot and the leads. In other words, 
we consider three different fluctuating in time potentials 
Vj{t), where the index j labels the left lead, the dot and 
the right lead. 

The wave functions are then modified in two 
ways. First, they acquire additional phase factors 

exp |i j^dt' eVj{t')^ = exp{(y9j(t)}. These phase factors 

can be eliminated by a gauge transformation, and we 
will omit them in what follows. In addition to that, 
the phase of the wave function changes every time the 
electron crosses one of the junctions. The corresponding 
additional phase factors acquired by the electron are: 



e~*¥'i, left junction, from left to right, 
e*'''^ left junction, from right to left, 
right junction, from left to right. 



where the matrix A reads 



right junction, from right to left. (A2) 



Below we will restrict ourselves to the case i > 0. This is 
sufficient, since the evolution operator at negative times 
is just a conjugate operator. Taking the phase factors 
(|A2I) into account we find U in the presence of fluctuating 

voltages Vj^{t) : 

U"^^ (t, 0; X, y) = A{t ~ x)U{t, 0; x, y)A*{-y) (A3) 



A{t) 











(A4) 



The relation (|A3|) turns out to be remarkably simple since 
every classical path - though being scattered arbitrarily 
many times inside the dot and from the outer barriers - 
can cross these barriers only twice: when it enters and 
leaves the dot. The times of these two events are related 
in a trivial way to the total time of the evolution and 
to the initial and final coordinates: For any x and y the 
path enters the dot at a time —y and leaves it at a time 
t — x. We also note that this property holds only for struc- 
tures with two barriers. In systems with three and more 
barriers no such a simple relation between If^ {t, 0; x, y) 
and U{t,0; x,y) exist. 



APPENDIX B: EFFECTIVE ACTION 

Let us now use the above expressions for the evolution 
operator and derive the effective action of a quantum dot. 
We start from the term Sr which is defined by the second 
term in Eq. Hll|l as 

iSR = ej dt' J rf3r(Gv+ai(t'-0,t',r,r) 

+ Gv+Mt' + 0, t\ r, r))V- [t\ r) . (Bl) 

Here the space integrals are taken over the dot and both 
leads. As we have already discussed, we will adopt the 
model assuming that the fluctuating voltage fields inside 
the leads do not depend on coordinates. Making use of 
Eqs. (I12I13II let us transform Eq. HB1|) to the following 
form 

iSn = 2i y dt' j d^rd^r'd^r"U^'' {t'Q,rr') 

X po(r', r")U^^ (Ot', r"r)^-{t', r). (B2) 

Let us split this expression by distinguishing each of the 
three coordinates r, r' and r" to be in the left lead, in the 
dot or in the right lead. Altogether we obtain 3'^ = 27 
terms. Without loss of generality this number can be 
reduced substantially by means of the following steps. 
First, we choose the initial density matrix po in the form 
corresponding to zero transmissions of both junctions. 
For such a choice initially we have three isolated systems, 
in which case the coordinates r' and r" always belong 
to the same electrode. This step already reduces the 
total number of terms down to nine. Second, we restrict 
ourselves to the limit of sufficiently long times, which is 
of a primary interest in our problem. In this limit the 
initial form of the density matrix inside the dot does not 
matter at all, electron transfer and relaxation processes 
will eventually lead to some final density matrix in the 
dot which will be independent of the initial one. Hence, 
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in the interesting limit of long times our effective action 
should not depend on the initial density matrix inside 
the dot and we can safely exclude the dot region from 
the integration over r' and r" extending this integration 
only at the two reservoirs. After that we are left with 
only six terms^, and the action Sr takes the form 



iSu — 2? 



E E f dt' [ d^r'd^r"pi{r'y') 
^ dPrUt!, it% rr')uf^ (Ot', r"r)^; {t'). (B3) 



In order to proceed further, we will make use of the 
unitarity of the evolution operator. It yields 

^ / d^rf/fj (t'Q, rr')uf^ {Qf , r"r) = 5{r' - r"!jB4) 



With the aid of this identity one can express the integrals 
over r in the dot region through the integrals over the 
coordinate in the leads, where an explicit form of the 
evolution operator l|A3p is already available. We find 

+ 2i f dt' J d^dy p'o{r',r") 

''d'rU^,Uot',r"r)Ur,Ut%rr')^-L{t') 



d^rUt,, (Oi',r"r)f/;,, (t'O, rr')^;,,(0 



where ipj^ — ip^ — ip~i, ~ Vd ^ V'r ■ The first term 
in this equation is proportional to the total number of 
electrons in the leads. Because of the charge neutrality 
this term should be exactly canceled by the analogous 
term from the ion background. For this reason we will 
omit this term in our further consideration. 

In order to evaluate the remaining terms we switch to 
the channel representation. Since we assume that the 
density matrix in the leads corresponds to the equilib- 
rium Fermi distribution of electrons, both this density 
matrix and the matrix of the phases ip~ are diagonal in 
the channel space 



Po{x -y) = Po{x - y)i, (p{t') 



i^R{t') 



The action Sr, acquires the following form: 



(B5) 



iSR = —2i / dxdytr 



Mx-y) / dt' / dzA{^y) 



X U{0,t';y,z)A*{t~ z)^-{t')A{t~ z) 



X U{t' ,Q]z,x)A*{-x) 



(B6) 



We also notice that the matrices ip{t') (jB5|) and A (|A4|) 
commute, hence A*(t — z) and A{t — z) drop out. Per- 
forming the change of the integration variables x —x, 
y — > —y, we obtain 



iSr — —2i / dxdy ti 



A*{x)po{y-x)A{y)B{t,y,x) 



Here we defined the matrix 



B{t,y,x) 



(B7) 



dt' J dz U{0, t'; -y, z)^- {t')U{t' , 0; z, -x). 

(B8) 

With the aid of Eq. (|Aip this matrix can be evaluated 
explicitly. For our purposes it is sufficient to set ip^ (t) = 
ip-{0) = 0. Then we get 



Bit,y,x) = 9ix)9iy) / dz [6iz ^ y)ip~ iz)Siz - x) 
Jo 

-S+{z~y)ip-{z)Siz-x)]. (B9) 



Making use of this expression together with (|A4I) and 
()B5p one can now directly multiply matrices in Eq. ljB7|) 
and arrive at the final result for Sr presented in Eq. H18|l . 

Let us now turn to the expression for Si which is de- 
fined by the last term in Eq. Hll|l 

Si = e2 f dt' f dt" /d3r'd3r"Gy+,i2(t',i",r',r")) 
Jo Jo J 

X V-{t", r")Gv+,2i(i", t', r", r'))V-{t', r'). 

The whole consideration is completely analogous to that 
carried out above for the term Sr. Making use of Eq. 
(O, we find 



Si 



e~ I dti j dt2 
10 Jo 



d^Xid^X2d^yid^y2d^Zid Z2 



X U^+{ti0,xiyi)po{yi,y2)U^+{0t2,y2X2}V (^2,2^2) 

X (^20, X2Z2)ho{z2,Zi)U^+ {Oti , Zi , Xi)V^ {ti , Xi) . 

where ho{z2,Zi) = 1 — po{z2,Zi). As before we identify 
the contributions containing the initial density matrix in 
the dot. Such terms can again be omitted in the inter- 
esting limit of long times. We again apply the unitarity 
condition (|B4I) and explicitly introduce the channel in- 
dices. Then after some straightforward manipulations 
we obtain 

Si — dti / dt2 / dxidx2dyidy2dzidz2 
Jo Jo J 

X tr|c/(0, ti;zi,xi)<f-{ti)U{ti,0; xi,yi)A* {-yi) 

xpo(yi - y2)A{-y2)U{0,t2;y2,x2)(p~ {t2) 

xU{t2,0; X2, Z2)A* {-Z2)ho{z2 - zi)i(-zi) \. 



18 



With the aid of the matrix (jB9|) this expression can be 
transformed further and eventually takes the form 

Si J dyidy2 Jdzidz2tT^^B{zi,yi)A*{yi)po{y2-yi) 

xA{y2)B(y2,Z2)A*(z2)ho{zi-Z2)A{zi)y (BIO) 

Multiplying matrices in Eqs. (jB10|l we arrive at the final 
expressions for Sj defined in Eq. (|21|l . 

APPENDIX C: AVERAGING OVER CIRCULAR 
ENSEMBLE OF 5" MATRICES 

Here we will evaluate the following average 

W = tr {ASHEi)BS{E2)CSHE3)DS{E4)). (CI) 

Let us use the representation (|38|l and split the 5— matrix 
into two parts 

S{E) + 5S{E), SS{E) = f'[l - U{E)]~^U{E)f . 

(C2) 

One can show that {5S{E)) — 0. At any given energy E 
the matrix tj (E) belongs to the circular ensemble. 

In a semi-classical limit averaging of the products of 
such matrices taken at the same energy can be performed 
with the aid of the diagram technique-^. Here we need a 
more general version of this technique which would also 
allow to average the products of matrices taken at differ- 
ent energies. Appropriate modifications of the rules can 
be formulated with the aid of the results^*. In short, 
using the terminology^^, every f/— cycle involving two 
matrices UiEi) and U'^{E2) carries an additional fac- 
tor M/{M — iEi2to) as compared to the situation^^. 
Here M = Nl + Nr, where A^l(_r,) is the number of 



channels in the left (right) barrier/lead, = 2tt/6 and 
Eij = Ei — Ej. An J7— cycle involving four matrices 
UHEi), U{E2), UHEs) and (/{Ei) gives an additional 
factor 

M^(M -i{E23+E4i)to) 

{M - iEiito){M - iE2ito)iM - iE2sh){M - iEi^to) ' 

Higher order J7— cycles do not contribute in the semi- 
classical limit M ^ oo, and we will not consider them 
here. 

Employing Eq. ljC2p one can express the average (|Cip 
as a sum of 10 non-vanishing terms. The first of them 
is trivial and equals ix{Aw BRCW DR). Several terms 
are bilinear in the matrices 5S and SS'^ . Such terms are 
evaluated with the aid of the formula 

M -trR'R' -iEi2to 

_ 2{trf'^Af){trfBf^) 

5(1 -lE.^rn) ^ ' ^ ' 
Here we have used the identity M—trR' R' — tvTT^ — g, 
and defined the dwell time td — 2to/g. Eq. I|C3|I is a 
direct generalization of the analogous formula^® derived 
for the case Ei — E2- This formula allows to evaluate 
the averages (EH). For example, = ti{P{E)i{E + 
cj)) = ti{CiS^E)C2S{E + w)), where we defined the 
matrices 

Employing Eq. (|(y3(l we arrive at H41|l . 

The terms involving three matrices SS and/or SS'' are 
slightly more complicated. They read 



f /4AQ^p ^RAQ^F ^r'AQtrP ^^ 4{trf'^Ar){tifcf^){trfBrR'^) 
tT{AdS(Ei)BSS{E2)CSS^{E3)) = ^ — — — ^ , 

g^[l - lEi3TD)[l - 1E23TD) 
g^{l - lEi2TDj{l - lEisTDj 



The derivation of Eqs. (|C5|) is straightforward, and we will not discuss the corresponding details here. 
What remains is to evaluate the average 

w = ti{ASS''{Ei)BdS{E2)CdS^{E3)DdS{Ei)). (C6) 

This average is given by the combination of five diagrams shown in Fig. El 

w = Di+D2 + D3 + Di-D5. (C7) 



The diagram 1)5 is subtracted in order to compensate for double counting of certain graphs in the diagrams D3 and 



19 



TAT* 



T'*DT' 



T'*DT' 




TAT 



T'*BT' 



TAT* 



TCI* 



T'*BT' 



T'*DT' 




TCT* 



T'*BT' 



TCT 



FIG. 3: Five diagrams contributing to tT{ASS^Ei)I35S{E2)CSS^E3)DSS{Ei)) . We refer to Ref . ^ for the definition of tfie 
graphical objects presented here. 



Di. The contributions of the diagrams read 



(trf^)(trFc)(trFBFz,) 
(M -tEiita){M ~iE23to)'' 
(trFAT(7ft)(trFB)(trFzj) 
{M - iE2ito){M ~ iE43to)' 



Da 



(trf^)(trfB)(trfc)(trFo)(M - ijE^s + ^41)^0) 
{M - tE^ita){M - iE2ito){M - iE2zh){M ~ iE^^to) ' 



{trfAf^Fc){tTFB){tYFD) „ (trfiftfCTt)(trFB)(tr/-D) 



{M-iE2ito)(M-iE43to)' 



{M ~ iE2ito){M ~ lEisto) 



Here we have introduced the ladder blocks defined in Fig. 01 



+ ftrff'fi" 1" / 

Fa = TiTt + y"(tTfAf'')R'R' \ ^ \ = tit^ + — 



5(1 - iEhtd)' 



Fh 



T' BT' + — -R'R', Fc^ TCT^ + — -R' R'\ Ed ^ T'' DT' + — -R'R'\ 

g(l-lE2lTD) g{l-lE23TD) g{l-lEi3TD) 
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Eventually we arrive at an important - though rather lengthy - formula 



2 (trf' ' i?f ') {tifCB^jDRAf^ ) 2 {trf'^BRCR''Df') (trf if t ) 



W 



g{l - iE2iTD) .g(l - iEiiTo) 

2{trf'^ DRAR'f Bf'){trfcf^) 2{tTf'''Df'){trfAR''BRCf'<) 



g{l - iE2?,TD) 



+ 



+ 



4(trf'^i?fO(trfC'ft)(tri?T'^i)i?ift) A{ivf'I3f'){irfAf'^){ivfCR'^bf'R') 
g''{l-iE2iTD){l-iE23TD) ^ g^{\-iE2iTD){l-iEiiTD) 

4(trf Z)f ') (trf if t ) {iiR'f'^BRCf^ ) 4(trf fif ') (trf Cf t ) (trf ii?t BTR'"^ ) 
g^{l-iEiiTD){l-iEi:iTD) ^ g^{l - iEi3TD){l - iE2^TD) 

4(trfift)(trf(7ft)tr(f'^i3f'f'^f)f') 4(trf '^Sf ')(trf '^flf ')(trf if tf Cf t) 
52(1 - iEiiTD){l ~ iE23TD) 3^(1 - iE2iTD)il - iE^sTo) 

8(trfift)(trfC'ft)(trf'^f)f')(trf'^i?f'i?'^i?') 8(trf if t)(trf Cf t)(trf '^Sf ')(trf '^flf 



g3(l - iE43TD){l - iE4iTD){l - iE23TD) 



.9^(1 - iE2lTDKl - lEiiTD){l - iE23TD) 



8(trf'^Bf')(trf'^f>f')(trfift)(trfC'ft7?'E'^) 8(trf '^i?f ')(trf '^flf ')(trf Cf t)(trf if +7?'^?'^) 



g^{l - iE2iTD){l - iEi3TD){l - iEnTo) g^{l - iE2iTD){l - iEi3TD){l - «i?23T£)) 

16(trf if t)(trfC'f^)(trf'^i?f')(trf'^f>f')(M - tvR'R'^R'R'^ - i{E23 + £^41)^0) 

g4(l - iE2lTD)il - i£^43TD)(l - «£^4lT£))(l - iE23TD) 



(C9) 



One can demonstrate that the result ljC9p is consistent 
with the unitarity of the S"— matrix. For this purpose 
let us put E4 = E3 and D — 1. Due to the unitarity of 
5'— matrices of the barriers almost all of the terms cancel 
out and the result reduces to 



W = tr AR^IBRC 



2(trf'^i?f')(trfC'ift) 



g{l - iE2iTD) 



= ii{AS\Ei)I3S{E2)C). 



(CIO) 



The last equation is a consequence of Eq. ljC3|) . This 
is just the expression expected from the definition 
and the unitarity of the 5— matrix. 

With the aid of Eq. ljC9|l one can find the averages v^j 
defined in Eq. H29(l . We transform them to the form 



= tv{CiS\E)CiS[E + uj)C2S\E + Lo)CiS{E)) 

= tr{CiS\E)C2S{E + uo)CiS\E + u)CiS{E)) (Cll) 



and similarly for and v^^^'. Making use of Eq. (jC9| 
we arrive at Eqs. H42() . Analogously, we derive the ex- 
pressions (|101|l . 



APPENDIX D: GENERAL EXPRESSIONS FOR 
THE CURRENT 



In Sec. 4A we focused on the expression for the current 
in the limit of zero external impedance Zs 0. Our 
analysis also allows to establish more general expressions 
valid for arbitrary Zs{uj). We obtain 

IiV)=IoiV) + SlLR + SlLL + 6lRR, (Dl) 

where 
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Io{V) = 2e f dxpo{x)T{x) (^e-Fr.L(x) ^^eV^x _ ^-F^,n{x) ^-^eVnx'^ ^ ^2) 



= 2elm J dx J dyidy2dzidz2 [pa{y2 - yi)ho{zi - Z2) + po{yi - y2)ho{z2 - zi)] 

X [(KLBi-yi) + KLL{-zi))i'\y2)f{z2) - {KRR{-y^) + KRL{-zi))?\y2)i{z2)\ } , (D3) 

(5/iL - -2elmjdx J dyidy2dzidz2 poiyi - y2)ho{z2 - zi)TLLiyuy2, z^, Z2)e'^''^'<-y'-y'-''+^^^^ 
X Tr{<t(a, + + - KLLhzi)){Siy2)S{z2)i - fHy2)f{z2)) 

+ (KRLi-yi) - KRLi-zi))i\y2)iiz2)]}, (D4) 

SIrr = 2elm J dx J dyidyzd^id^s Po(yi - y2)/io(z2 - 2i)^i?^fl(2/i, ?/2, ^^i, Z2)e-*^«(^i-^^-"^+"^)l 

X Tr{f'\a: + zi)i'ix + yi) [{KRR{-yi) - KRR{-zi)){5{y2)5{z2)i - r'\y2)f' {7.2)) 

+ {KLRi-yi) - KLR{-zi))i'\y2)i'{z2)] }. (D5) 

Here we defined 

^LF!{yi,y2, Zi, Z2) = c^^rLrL^y^^y^)^^i-i-^^-^^^-'^)+^LR{yi-zi)+FLR{y2-z2)-FLR(yi-z2)-FLR(y2-zi) ^ (J^&) 

J^LL{yi,y2, zi, Z2) — c^^^'^^y^~y^)~^^^'^'^^^'''^'>^^^'^'^y^~''^)^^i^i^^y^-'^^)+^^'^'^vi~z2)+FLL(y2-zi) ^ ^y)T) 

/duj CO 
— ImK.j (lu) coth — (1 - cos 
ZTT 21 

^ "''^ \ ± / k,n=L.,R 

5vkn{'-^) are defined in Eq. I|29l) . The function Trr is obtai ned by substituting Fll{x) Frr{x) into Eg. I|D7|). 
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